Structural, orbital, and magnetic order in vanadium spinels 
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Vanadium spinels (ZnV204, MgV^CU, and CdV204) exhibit a sequence of structural and magnetic 
phase transitions, reflecting the interplay of lattice, orbital, and spin degrees of freedom. We offer 
a theoretical model taking into account the relativistic spin-orbit interaction, collective Jahn- Teller 
effect, and spin frustration. Below the structural transition, vanadium ions exhibit ferroorbital 
order and the magnet is best viewed as two sets of antiferromagnetic chains with a single-ion Ising 
anisotropy. Magnetic order, parametrized by two Ising variables, appears at a tetracritical point. 



Frustrated magnetic systems have generated interest 
among researchers for more than a decade Q, 0, Q ■ Frus- 
tration refers to the inability of spins to satisfy conflicting 
interactions resulting in a large degeneracy of the ground 
state at the classical level. The effect is particularly 
strong for spins with antiferromagnetic interactions on 
the three-dimensional network of corner-sharing tetrahc- 
dra, known as the pyrochlore lattice |4| . Magnets of this 
kind have rather unusual properties. For example, spinels 
containing magnetic ions Cr 3+ (ZnCr 2 04, MgCr 2 C>4, and 
CdC^CU) undergo a magnetoelastic phase transition 
0- , even though they have no spin chains commonly 
associated with the spin-Peierls effect 0. 

In this Letter I discuss magnetic and structural prop- 
erties of another remarkable group of magnetic spinels 
containing V 3+ ions (ZnV 2 4 , MgV 2 4 , and CdV 2 4 ). 
These compounds represent a step up in complexity from 
their chromium counterparts because vanadium ions have 
an orbital degree of freedom — in addition to lattice vi- 
brations and spin. The interplay of orbital, spin, and 
vibrational motion against the backdrop of a high spa- 
tial symmetry leads to a sequence of phase transitions 
involving changes in structural and magnetic properties. 
The physics involves the relativistic spin-orbit interac- 
tion, collective Jahn- Teller effect, and spin frustration. 

The vanadium spinels undergo a structural phase tran- 
sition at a temperature T s and a magnetic ordering at a 
lower temperature T/y. The structural distortion low- 
ers the crystal symmetry from cubic Fd3m to tetrago- 
nal IAi/amd HQ, with a slightly shorter lattice period 
along the c axis. The Neel temperature Tjy is rather 
small compared to the Curie- Weiss scale 0cw, a sure sign 
of strong spin frustration. Ordered magnetic moments 
point along the c axis and have the length (/i) < fig at 
helium temperatures. Magnetic susceptibility is history- 
dependent both below and above T/v [jj, 0] . The num- 
bers for ZnV 2 4 are re pres entative of the entire family: 
c/a = 0.994, T s = 51 K 10] , T N = 42 K H2, 9 C w = 420 
K [13, and (n) = Om^Tm- 

The main difference between Cr 3+ and V 3+ ions is the 
number of electrons in the partially filled 3d shell. The 
crystal field splits the 3c? quintet into a low-energy triplet 
t2 g and a high-energy doublet e g , separated by a gap of 
roughly 2 eV. In accordance with Hund's rules, the three 



FIG. 1: Orbital order of vanadium ions below the structural 
phase transition. A view down the [001] direction. Smaller 
symbols are farther away, (a) The result of Tsunetsugu and 
Motome 0. (b) This work. 



d electrons of Cr 3+ have spin 5=3/2 and occupy all 
three triplet levels. The ion has a spin degree of freedom 
only. In contrast, the two d electrons of V 3+ have spin 
5 = 1 and occupy two out of three ti g orbitals. 

Tsunetsugu and Motome have recently developed 
a model explaining structural and magnetic properties of 
the vanadium spinels. Following the approach of Kugel 
and Khomskii, they focus on Coulomb and exchange in- 
teractions between magnetic ions. However, their model 
is at odds with experimental observations: the orbital 
order predicted by the theory is incompatible with the 
spatial symmetry of the tetragonal phase, identified as 
IAi/amd || 0- The staggering of ions with empty yz 
and zx orbitals would break reflections m in the planes 
(110) and (110) and diamond glides d in the planes (100) 
and (010), lowering the symmetry to I4i/a [Fig.^a)]. 

I propose an alternative theory that is compatible with 
the structural data and predicts rather nontrivial mag- 
netic behavior. In the new model, the structural order 
arises from simple single-ion physics: the interplay of the 
JT effect and relativistic spin-orbit interaction. The i 2g 
orbitals are partially ordered: their occupation numbers 
satisfy n yz = n zx < n xy [Fig. IHb)] . The distorted mag- 
net is best decribed as a collection of antiferromagnetic 
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chains, running in the directions [110] and [110], with an 
effective spin J = 2 and a single-ion anisotropy of the 
Ising type. Parallel chains are weakly coupled; the cou- 
pling of orthogonal chains is frustrated by symmetry. As 
a result, the Neel phase is characterized by two indepen- 
dent Ising order parameters; the magnetic transition oc- 
curs at a tetracritical point. Two nearby magnetic phases 
exist in which only one half of the chains order. Finally, 
memory effects in magnetic susceptibility are tied to the 
spontaneous breaking of the lattice symmetry. 

Choice of the model. The main difficulty in obtaining 
a consistent theoretical picture of the vanadium spinels is 
a poor separation of energy scales. The ti g orbitals of a 
B site in spinels point away from the oxygen ligands and 
thus have a relatively weak JT coupling [l5j comparable 
in magnitude to the relativistic spin-orbit coupling. (A 
recent estimate 0] puts the JT energy in YVO3 at 20-60 
meV. The spin-orbit splitting of the 3 Fi and 3 J 7 3 levels 
in a free V 3+ equals 39 meV [13 •) Furthermore, the 
effective interaction of spin and orbital degrees of freedom 
on adjacent vanadium ions has a similar energy scale [l3j ]. 
It is not at all obvious which of the three interactions 
plays the dominant role. 

To make progress, one can adopt a phenomenological 
approach and investigate three simple limits, in which 
the dominant interactions are, respectively, (I) the JT 
coupling, (II) the spin-orbit cou plin g, (III) the Kugel- 
Khomskii interaction. Model III 13] yields orbital order 
incompatible with the observed symmetry of the tetrag- 
onal phase H, E| . We therefore explore limits I and II. 

Basis states and notation. If the energy difference be- 
tween the 3 P and 3 F electronic terms of the 3d 2 con- 
figuration (w 1 eV) is neglected in comparison with its 
crystal- field splitting in a VOg octahedron (w 2 eV) , the 
3 Ti g ground state of a V 3+ ion has the t\ g e® g configura- 
tion. Its orbital part is triply degenerate and has the 
following basis: 

\X) = (\zx)\xy)-\xy)\zx))/V2, 

\Y) = (\xy)\yz)-\yz)\xy))/V2, (1) 

\Z) = (\yz)\zx)-\zx)\yz))/V2. 

In this Hilbert space, the operator of orbital momentum 
L can be expressed in terms of an effective angular mo- 
mentum L' of length L' = 1 such that 

L' z \X)=i\Y), L' z \Y) = -i\X), L' z \Z)=0 (2) 

etc. It can be verified that L = —ctL', where a = 1 
in the approximation made Inclusion of electronic 
correlations mixes in the states t^e* and t%„e 2 an d puts 
a somewhere in the range between 1 and 3/2 "Jig. 

Model I. A tetragonal distortion of a VC>6 octahedron 
along the c axis splits the 3 Ti ff state into a doublet and 
a singlet. The splitting can be mimicked by the traceless 
operator 

V JT = 8[L' 2 z -L'(L' + l)/3], (3) 



where 8 is the tetragonal strain expressed in energy units. 
For a 3c? 2 ion, 8 > for an elongated octahedron [l5| . 
The lowest-energy state is \Z), in which the xy orbital is 
empty and the energy is —26/3. In the opposite case (a 
flattened octahedron, 8 < 0), the orbital ground state is 
doubly degenerate and has a higher energy — 1<5|/3. Inclu- 
sion of the elastic energy 0{8 2 ) and minimization with 
respect to 8 yields a ground state with 8 > 0. Thus Model 
I, dominated by the JT coupling, predicts an elongation 
of V0 6 tetrahedra contrary to the data [1 H Eol [Til Hof . 
Model II. The relativistic spin-orbit coupling in V 3+ 

V LS = A(L • S) = -a\(L' ■ S) (4) 

has the strength A ~ +20 meV [l^- The "total angular 
momentum" J' = L' + S is a conserved quantity. The 
lowest-energy state is the J' = 2 quintuplet separated by 
the energy gap 2aA from the J' = 1 levels. 

Inclusion of the JT interaction J3J as a perturbation 
(8 <C A) splits the J' = 2 quintuplet into three levels. 
With the aid of the Wigner-Eckart theorem, one finds 

V 3T = {5/2)[J'l-J\J' + l)/3]. (5) 

Remarkably, the ground state has the same energy — 1<5| 
regardless of the sign of the 8. In an elongated VC>6 oc- 
tahedron (8 > 0), the ground state is nondegenerate and 
has J z = 0; its magnetic moment — hb{— «L' + 2S) is 
quenched. In a flattened octahedron (8 < 0), the ground 
state is a non-Kramers doublet J' z = L' z +S z = ±2. Inclu- 
sion of the elastic energy 0(S 2 ) fails to determine the sign 
of the tetragonal JT distortion. It will be determined by 
several competing perturbations, such as higher-order JT 
terms, Coulomb and spin-exchange interactions between 
vanadium ions. 

The effect of Coulomb interaction can be crudely esti- 
mated by counting the number of filled ti g orbitals with 
the largest overlap. For example, the Coulomb energy of 
a [110] V-V bond is raised if the xy orbitals are occupied 
in both ions |l3l |. In a state with uniform orbital order, 
the average occupations of the orbitals are as follows: 

J' z = ±2: (n yz ) = (n zx ) = 1/2, (n xy ) = 1, , . 

J' z = : (n yz ) = (n zx ) = 5/6, (n xy ) = 1/3. W 

The expectation value of the Coulomb energy, summed 
over nearest-neighbor V-V pairs, will be proportional to 
(n yz } 2 + (n zx ) 2 + (n xy ) 2 . This expression evaluates to 3/2 
for both J' z = and J' z = ±2 ground states. Thus the 
Coulomb repulsion between V 3+ ions is the same in the 
elongated and flattened ground states of VOq octahedra 
and docs not influence the outcome of their competition. 

(The accidental degeneracy of the Coulomb term can 
be traced to the degeneracy of the flattened and elon- 
gated JT states for J' = 2. Both the Coulomb and 
JT energies can be expressed in terms of the operators 
J'l - J' {J 1 + l)/3 and J'l - J' 2 It is therefore likely 



that the conclusion reached in the previous paragraph is 
not sensitive to the approximations made in this paper.) 

Exchange energy, on the other hand, clearly favors the 
magnetic ground states J' z = ±2 because it can be made 
negative by an appropriate choice of spin orientations. 
In contrast, exchange energy is zero if the ions are in the 
nonmagnetic state J' z = 0. Anharmonic terms of the JT 
energy 0(5 3 ) can favor either sign of S. Therefore, one 
cannot deduce the sign of the distortion on general the- 
oretical grounds. Nonetheless, experimental observation 
of antiferromagnetic order in ZnV2C>4 at low tempera- 
tures points to the magnetically active states J' z = ±2 
characterized by a flattening of the VC>6 octahedra along 
the c axis. For a given direction of the spin, S z = ±1, 
one of the electrons occupies the xy orbital, whereas the 
other is in the state \yz) ± i\zx), whose actual orbital 
moment L z = —aL' z = =fa points opposite to the spin. 
Equal occupation of the yz and zx orbitals in such a 
state is in agreement with the observed spatial symmetry 
of the crystal (IAi/amd). Magnetic moments of length 
(2 — a)fJ-B, where 1 < a < 3/2, point along the c axis, as 
found experimentally. Thus phenomenological Model II, 
with strong relativistic spin-orbit coupling, a moderate 
JT distortion, and weak V-V interactions, is compatible 
with experimental observations. 

Structural order and symmetry-breaking fields. Below 
T s the VC>6 octahedra flatten along one of the three ma- 
jor axes. From the symmetry viewpoint, the long-range 
order is similar to that of the 3-state Potts ferromagnet. 
A two-component order parameter f = (/i,/a) can be 
constructed, e.g. from the lattice constants: 

h = (a + b - 2c)/ V6, h = {a- b)/V2. (7a) 

Other observables can be used to define the order param- 
eter. The lattice constants (a, 6, c) can be replaced with 
orbital populations ((n yz ), (n zx ), (n xy )), or with mag- 
netic susceptibilities Xij — dMj/dHi: 

fl = (Xxx + Xyy-2Xzz)/V6, h = (Xxx~ Xyy) I 'V2. (7b) 

Explicit construction of the order parameter allows us 
to identify symmetry-breaking fields that can be applied 
to select a particular phase with broken symmetry. The 
most obvious among these is the two-component (tetrag- 
onal and orthorombic) stress h = (hi, h^)'- 

hi = (e xx + e yy - 2e zz )j V6, hi = (e xx - e yy )/V2. (8a) 

Cooling under a weak uniaxial stress applied along the c 
axis will produce an ordered phase with a = b > c. 

A less obvious symmetry-breaking perturbation is 
magnetic field H. The magnetic energy density of a sam- 
ple can be written as -(xxxH x +x yy Hy+XzzH z )/2. Un- 
less the field is parallel to a (111) direction, it lowers the 
energies of different ordered phases by different amounts. 
The symmetry-breaking field is 

hi = {Hi + Hi - 2H*)/V6, h 2 = (H 2 X - H*)/V2. (8b) 
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FIG. 2: Two of the four magnetically ordered states. The 
signs refer to the two Ising order parameters (o-[ 110 j , crnjo])- 
Ground states (a) and (b) are related by a (110) reflection m 
or by a (110) glide a in a (001) plane. 

These considerations could explain memory effects in 
magnetic susceptibility of ZnV2C>4 above the Neel tem- 
perature Tn 0, E3 ■ A crystal cooled through the struc- 
tural transition in magnetic field chooses the distorted 
phase with the lowest magnetic energy. The difference 
between %fc and Xzfc in powder samples may simply 
reflect the presence of the order parameter l|7b|) below T s 
and its sensitivity to the symmetry-breaking field (|8b|) . It 
would be highly desirable to study these memory effects 
in single crystals. 

Anijerromagnetic Ising chains with frustrated coupling. 
In the distorted state, the magnetic system, comprising 
spins and orbital moments, can be thought of as a col- 
lection of chains running in the [110] and [110] direc- 
tions (Fig. I2J) . Because the xy orbitals are always occu- 
pied, the largest exchange coupling takes place along the 
chains. Exchange coupling along the other (110) direc- 
tions is expected to be substantially weaker because the 
relevant orbitals (yz and zx) are only half-filled. Onoda 
and Hasegawa Q have extracted the values of J[no] ~ 5 
meV and J[on] = J[ioi\ ~ 3 meV for CdV2C>4. In ad- 
dition to the single-ion Ising anisotropy J5j, which is a 
product of orbital order and spin-orbit coupling, the ex- 
change interactions are likely anisotropic as well. 

As the distorted magnet is cooled further, spin corre- 
lations develop first along the chains, as indeed observed 
[lflj . If the interchain coupling is weak, long-range mag- 
netic order will set in at a low temperature, when each 
chain already has a large correlation length. In that case, 
the magnetic transition will have almost no effect on the 
energetics of the spins. Indeed, measurements of specific 
heat in MgV2C>4 reveal a very weak singularity at the 
magnetic transition (1 1| . 

While magnetic interactions between parallel chains 
are expected to be weak because they are well sepa- 
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FIG. 3: The temperature-strain phase diagram in the vicinity 
of the tetracritical point T = Tjv, e xy = 0. Gray symbols 
correspond to disordered spins. Thick solid lines are phase 
boundaries. Other notations are same as in Fig. [2] 



tions, orbital motion, and spin. It is desirable to put 
these considerations on a quantitative footing with the 
aid of a microscopic theory. Nonetheless, the qualitative 
predictions of this work — a spin gap due to single-ion 
anisotropy, a magnetic tetracritical point, and a poten- 
tial relation between magnetic memory effects and broken 
lattice symmetry — can be tested experimentally. 

Upon completion of this work I learned that a similar 
model had been considered by D.I. Khomskii. 
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National Science Foundation under the Grant No. DMR- 
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rated in space, the coupling between crossing chains may 
be substantial: four of their spins are nearest neighbors 
forming a tetrahedron (Fig.|2J|. However, in this case the 
coupling is frustrated by the lattice symmetry: a (110) 
reflection m (or a glide a) changes the sign of staggered 
magnetization on the chain running in the direction [110] 
but preserves staggered magnetization on the [110] chain. 
Therefore, even if the crossing-chain coupling is numer- 
ically comparable to the intrachain coupling, it is still 
rendered weak by geometrical frustration. 

Magnetic order. Neutron scattering indicates that ad- 
jacent parallel chains have equal values of staggered mag- 
netization (Fig. |2J). Because of the geometrical frustra- 
tion, the system has four ground states: those shown 
in Fig. |21 and their time-reversed copies. These ground 
states can be characterized by two independent Ising or- 
der parameters: one staggered magnetization for each set 
of chains, <7[no] and cruigt [iol ]. 

Magnetic tetracritical point. Both sets of chains ap- 
pear to be ordered at the lowest temperatures [j| . In the 
simplest scenario, both Ising order parameters develop 
continuously at the magnetic transition temperature TV, 
which is characteristic of a tetracritical point [2lL l22| . 
One should be able to observe two additional phases in 
the vicinity of the tetracritical point, in which only one 
set of chains is ordered, while the other is not. This may 
be achieved by creating an orthorombic strain e xy in the 
vicinity of the magnetic transition. Doing so will violate 
the equivalence of the [110] and [110] chains while still 
preserving the geometrical frustration. A temperature 
scan under a uniaxial [110] stress should reveal two sep- 
arate magnetic transitions for each set of chains. The 
temperature-strain phase diagram in the vicinity of the 
tetracritical point is shown in Fig. |2| 

Outlook. A phenomenological approach taken in this 
paper is justified by the complexity of vanadium spinels: 
the relevant degrees of freedom include lattice vibra- 
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